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D 1 Abstract. We give a status report of convective Cepheid and 

C/3 , RR Lyrae model pulsations. Some striking successes can be 

reported, despite the use of a rather simple treatment of tur- 
bulent convection with a ID time-dependent diffusion equation 
for the turbulent energy. It is now possible to obtain stable 
double-mode (beat) pulsations in both Cepheid and RR Lyrae 
| models with astrophysical parameters, i.e. periods and am- 

. plitude ratios, that are in agreement with observations. The 

ON I turbulent convective models, however, have difficulties giving 

£^ . global agreement with the observations. In particular, the Mag- 

\ ellanic Cloud Cepheids, that have been observed in connection 

with the microlensing projects have imposed novel observational 
constraints because of the low metallicity of the MCs. 

6 

C3 , 1. Introduction 

Cepheid and RR Lyrae variables have played a central role in astrophysics. This 
^ ■ is not only because of interest in the pulsation mechanism itself, and of the diffi- 

culty one had to reconcile them with stellar evolution and stellar pulsation (Cox 
1980), but also, in the case of the Cepheids, because of the essential role they 
have played as a local distance indicator in the quest for determining the Hubble 
constant. As a result they are the best known and studied stellar pulsators (for 
a recent review of stellar pulsations, cf. Gautschy & Saio 1995). 

Recently, our observational base of Cepheids has been greatly enlarged with 
the data from the microlensing projects EROS, MACHO and OGLE (e.g. Fer- 
let et al.1996). Not only has the number of well observed Cepheids increased 
dramatically through these efforts, we now also have a much broader data base, 
particularly because the Magellanic Cloud (MC) Cepheids have metallicities that 
range from one half to as low as one quarter solar. 

This increased observational knowledge of Cepheids has also brought to 
light new problems in the modelling (e.g. Buchler 1998). For example, purely 
radiative models give light and radial velocity curves in good agreement with the 
observations of the Galactic Cepheids (Moskalik et al. 1992). But is no longer 
true for the Magellanic Cloud Cepheids with their low metallicities. 
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The reader may ask here why theorists have bothered to study purely ra- 
diative models (in which heat transport is assumed to be through radiation 
transfer only) when we know that convection plays a role in these stars, notably 
in diminishing vibrational driving in the colder models and in bringing about 
the red edge of the instability strip. The reason is that radiative models have 
been expected to give excellent light curves and radial velocity curves because 
convective transport was deemed to be not very efficient nor presumably impor- 
tant overall. And, indeed, this expectation was satisfied for Galactic Cepheid 
models. However, despite great efforts, it has not been possible to model the 
MC Cepheids with purely radiative models, and it is fair to say that all possible 
reasons for the discrepancies have now been exhausted (Buchler 1998). For a 
similar conclusion about RR Lyrae see Kovacs & Kanbur (1997). 

We have therefore been led to consider convective transport in our models. 
The physical conditions are such that one expects well developed turbulence: 
the Rayleigh number in the stellar envelope is huge and the Prandtl number is 
tiny. Possibly plumes play an important role (Rieutord & Zahn 1995, Zahn, this 
Volume) . 

We are not merely interested in computing a convective flux in a static 
model, we also want the linear eigenvalues in order to gain knowledge about 
the periods of the models and their stability. Finally, we want to be able to 
do hydrodynamics so as to compute the longterm pulsational behavior of the 
models, whether period or multiperiodic (as for the Beat Cepheids, q.v. below). 
Since we are interested in the pulsational problem which is complicated enough 
in itself we need a simple, but physically and mathematically robust recipe for 
the description of turbulence and convection. 



2. The Equations 

The Cepheids are radial pulsators in which the centrally condensed core is in- 
ert, and only the envelope undergoes pulsational motion. The motions are thus 
governed by ID hydrodynamics in spherical geometry, supplemented by a de- 
scription for the coupling of turbulence and convection with the fluid dynamics. 

Hydrodynamic Equations 



du 
~dt 

de dv 

hp — 

dt y dt 



Id, 

— Q~ (P + Pt+Pu ) 

p or 
1 d 



GM r 



pr 2 dr 



r 2 {F r + F c ) -C 



(1) 

(2) 



Convection interacts with the hydrodynamics of the radial motion through 
the convective flux F c term, through a viscous eddy pressure p u and a turbulent 
pressure pt, and finally, through an energy coupling term C. The question is how 
to approximate these quantities in the simplest physically acceptable way, so as 
to make it still possible to compute nonlinear stellar pulsations. 
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The simplest recipe involves a single, time-dependent diffusion equation for 
the turbulent energy e<, of the form 



Turbulent Energy Equation 



de t dv Id 

-Tr + {Pt + Pu) -jr = 

at at pr l or 



r 2 F t ) +C 



(3) 



The coupling term has the general form 

C = S-e 

In the spirit of Kolmogorov the dissipation rate e is taken to be 

e = a d e'l /2 /A 



(4) 



(5) 



where A = a\H p is the mixing length, proportional to the pressure scale height 
H p = d Inp/dr = p/(pg). We define the turbulent pressure pt = a p pet and the 
viscous eddy pressure p u = 4/3a v a\ H p ^felr{d{u/ r) / dr) . 

There is however no unique way of defining F c and S in terms of e t and Y, 
and several variants have been used. Stellingwerf (1982), Kuhfu6(1986), Geh- 
meyr (1992), Gehmeyr & Winkler (1992), Feuchtinger (1998a), Bono & Stelling- 
werf (1994), Bono et al. (1997, 1999) and Yecko et al. (1998, hereafter YKB) 
all have used such a ID equation in the computation of nonlinear pulsations, 
but they have made different possible choices for the dependences of the connec- 
tive flux F c and for the coupling term C on the turbulent energy et and on the 
dimensionless entropy gradient 



Y - Hp — 

c„ dr 



With the definitions 



A = a c a\ pCpT 
B 



p(3T (3T 
a s a\ \ = a s a\\ — c s 

VP V r i 



(6) 

(7) 
(8) 



where (3 = (dlnv/ dT) p and c s is the sound speed, we can write the three schemes 
as 



F c = Ae\ /2 Y 

S = a d B 2 e\ ,2 Y/K (GW) 
F c = Ae\ /2 Y 

S = a d Be t Y 1/2 /A (YKB) 
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F c = A/B e t Y 1/2 
S = a d Be t Y^/K 



(9) 



The combinations of A [energy /volume] (essentially the internal energy) 
and B [velocity] (essentially the sound speed c s ) have been chosen so that in the 
absence of diffusion all three schemes give the same static equilibrium model. 
(Generally they do not give same linear vibrational eigenvalues though). The 
model involves seven dimensionless a parameters of order unity. 

All three recipes can be parametrized with a coefficient N such that for 

N > : F c = A e\ /2 Y 

N<0: F c = A/Be t Y 1/2 (10) 

and, with b = 1/\N\, 

S = a d a s B 2b e 3 t /2 - b Y b /A (11) 

Thus we have Gehmeyr- Winkler (or Kuhfufi) : N=l, YKB: N=2, Stellingwerf: 
N=-2. (Note that one more unused possibilities exist, viz. N=-l). 

It is of course interesting to see the influence of the chosen recipe on the 
stellar pulsation, which we address in the next section. 

At this point we should note that of course more complicated recipes have 
been suggested, none of which have been implemented in nonlinear stellar pul- 
sations though. For example, in a much quoted, but unpublished 1968 preprint 
Castor reduced the problem of turbulent convection to a set of 3 coupled time- 
dependent diffusion equations for the three second order moments of the vertical 
velocity fluctuations w and the temperature fluctuations 6, namely < ww>, tur- 
bulent energy, < w6 > and < 66 >. The set of equations was closed with a 
down-gradient approximation for the respective fluxes, and with an expression 
for the turbulent energy dissipation e oc e\ /2 (Eq. ||). Kuhfu8(1986) also consid- 
ers a 3-equation version of the ID recipe. 

Recently, in a series of papers, Canuto (e.g. Canuto 1998) reexamined this 
problem and extended this formalism to higher order in which the downgradi- 
ent approximations are avoided because they have been found lacking, in which 
vertical-horizontal anisotropy is allowed for through the introduction of a sep- 
arate vertical turbulent energy, and finally, in which a dynamic equation for 
the energy dissipation e is introduced. As a result there are now five nonlin- 
ear coupled time-dependent diffusion equations to be solved together with the 
hydrodynamics and radiation transport equations. For the time being, incorpo- 
rating such a scheme into our pulsation code is a daunting task and we prefer 
to get as much insight from a simplified one-equation treatment as possible. 

It is nevertheless of interest to see to which of the above recipes Canuto's 
formalism reduces when additional, 'usual' simplifying assumptions are made. 
Thus, starting with Canuto & Dubikov's (1998, hereafter CD) Eqs. (CD19a- 
19d) we keep the first for K (= et) which is the equivalent of our equation for 
the turbulent energy after we make the down gradient approximation Df(et) oc 
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eJ /2 Vet. Instead of Eq. (CD19d) one can introduce an equation for the anisotropy 
turbulent energy by subtracting 1/3 x Eq. (CD19a) from (CD19d). In this equa- 
tion and in Eqs. (CD19b-19c) one now makes a local time-independent approx- 
imation, viz. d/dt — ► 0, and ignores the diffusion terms Dt and the x terms 
(large Peclet number). After some algebra this leads to 



F c = dipepT (e 2 /e) Y /V = d x pc p T e\ ,2 Y /V (12) 
S = dx (at/Tx) (e 2 /e) /V = d x (a t /T x ) c 2 s e] /2 Y /V (13) 

V =l-d 2 c 2 s (e t /e) 2 = 1 - d 2 a\ (c 2 s /e t ) Y (14) 

In each of these last three equations, the second expression is obtained with the 
use of Eq. The d\ and d 2 are positive dimensionless constants of order unity. 

Except for the appearance of a denominator T>, the CD formalism, with the 
additional quoted approximations, thus leads to the Gehmeyr- Winkler recipe. A 
quick inspection shows that the GW recipe is equivalent to ignoring anisotropy, 
i.e. to setting <w 2 >= | et, instead of using Eq. (CD19d), and to ignoring the 
<99> in Eq. (CD19c). 

This denominator, however, has a pole because d 2 is positive and, if used 
in a pulsation code, it would be disastrous. Indeed, for small et, Eq. ^ implies 
that below a certain threshold any existing et will decay away! The reason for 
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Figure 2. Integration of det/dt = C(et) with small initial et (bottom 
curves) and large initial ej (top curves) for the two values of \N\ in Eq. [ll] 

this unphysical denominator is either in the starting formalism, or more likely, 
in the additional approximations that we have introduced, 

The disregarded time derivatives of <98> and of J =<w8> (~ F c ) is not 
justified when ej very small, and can dominate when turbulence is about to set 
in. To see this, we look at the equations for small et when we need to consider 
only the source terms in Eqs. (CD19) (the dissipation terms having a higher 
power of et) which leads to e.g. 

c 2 

d 2 X/dt 2 = ga (3X = (aT/T^-^Y X = l/r 2 r X (15) 

p 

where, to within factors of order unity, X is any of the second order moments, 
and T gr is the initial growth-time of X. With this in mind we can write 
(l/X)(dX/dt) = (I /e t ) (det/dt) = (l/e t )gaJ (CD Eq. 19a). It is easy to see 
that when the time derivative dominates (at small et) then, <96 > instead of 
being proportional to J, one becomes 

<69>= 2(3(ga)K 

which thus removes the troublesome denominator. 

It is possible to carry out the intended plan of only keeping Eq. (CD 19a) 
and of replacing the three equations (CD 19b-19d) by their local limit. However, 
as we have shown it is necessary to keep the derivative terms at the cost of a 
higher order polynomial equation for J (F c ) which prevents simple analytical 
expressions such as Eq. ([13), and is not very practical. 

3. Comparison of Recipes 
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Figure 3. Light curves (L/L Q ) and radial velocity (km/s) curves for N=l 
and ±2 
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Figure 4. Variations of the turbulent diffusion and turbulent decay times 
(s) over the Cepheid envelope. 



3.1. Time-Dependence 

To give an indication how important time-dependence is in the pulsating Cepheid 
envelope we show the behavior of the turbulent energy as a function of time 
over a period for a typical cepheid model in Figure |l|. The region of maximum 
turbulent energy clearly undergoes substantial motion with respect to the fluid. 
It also shows a widening and intensification of the turbulent region during the 
pulsational compression phase. 

Because of this time-dependent behavior of the turbulent energy and be- 
cause of the differences in the expressions for the source terms in Eqs. 1C it is of 
interest to see how et responds in the various formulations. Gehmeyr and Win- 
kler (1992) and Kuhfuss (1986) already compared some of the properties of their 
scheme to that of Stellingwerf. Here we now show some additional comparisons 
of the various schemes. In Figure § we display the temporal responses of et to a 
constant entropy gradient Y . Thus we show the results of the integration of 



de t 
dt 



C 



(16) 



for two initial values of et, a very small one, and a large one. The temporal 
response is quite different for the two values of \N\, the adjustment time for 
et to a changing source Y being much slower for the larger N, i.e. for the 
Stellingwerf choice of source. A linearization of Eq. [l6| shows that the kernel has 
a different et dependence which can lead to small differences in the spectrum of 
the eigenmodes. 

3.2. Limit Cycles 

Because the different schemes have different timescales for the turbulent energy, 
it is of interest to compare the nonlinear behavior of models. In Fig. |3| we 
show the fundamental limit cycles (light curves and radial velocity curves) for a 
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Cepheid model with M=5M , L=1664L Q , T e// =5400K, X=0.70, Z=0.02, (a d 
= 4.0, a c = 3.0, «a = 0.35, a s = 0.75, a u = 1.5, a* = 1, a p = 0.67). Results 
are displayed for the three formulations N=l and N = ±2. The calculations 
actually show a great insensitivity to N, despite the expectations of the previous 
section. The reason for this insensitivity is that there are several time-scales at 
work, viz. a source time, a decay time and a diffusion time for ej. In addition, 
for each of the time-scales varies by many orders over the stellar envelope, as 
can be seen in Figure [| which shows the run of the diffusion and of the decay 
time-scales for the same Cepheid model. The insensitivity is also in accord with 
the results of YKB who showed that a 'sudden' approximation (C = 0), i.e. an 
instantaneous adjustment of &t to the local static equilibrium value gives quite 
a good approximation to the exact linear eigenvalues. 

3.3. Work Integrand 

It is well known that the pulsations of the classical variable stars are self-excited 
through the n mechanism (e.g. Cox 1980, Gautschy & Saio 1995). When, for a 
given vibrational mode, the driving that occurs in the partial ionization regions 
overwhelms the damping that occurs elsewhere, the mode becomes unstable 
or self-excited. This can best be demonstrated by the work-integrands which 
show the cycle averages of < pdv > throughout the star, i.e. the amount of 
internal energy converted into mechanical (pulsational) energy. The pressure 
that appears here is the total pressure p = p g + r + Pt + Pv We can thus look at 
the contributions of these components separately. 

In Figure [| the total work-integrand for linear perturbations is displayed as 
a thick line for the fundamental Cepheid model of the previous section. Shown 
with arrows are the points of half-ionization of H at 11,000 K and and He 7 
41,000 K, respectively. The molar tooth-shaped feature on the right is produced 
by the combined H and first He ionization regions which cause a broad convective 
region, with the spikes occurring at its boundaries. Just as in the purely radiative 
models the second ionization of He produces both a driving and a damping 
region. 

The turbulent pressure plays a very small role in general, although turbulent 
pressure gradients can be important locally (YKB). The eddy pressure, on the 
contrary, is an important source of damping, especially in the limit cycle where 
turbulent energy is spread over a wider range. For reference, the bottom figure 
shows the source of turbulence, namely the dimensionless entropy gradient Y, 
the turbulent energy, normalized to unity, and the fraction of energy carried by 
the convective flux F c . The driving is still positive in the H/He partial ionization 
region even though F c is large there. Some models, especially those with Z=0.02 
can also have a convective region associated with Fe-group elements. 

As is typical in the partial H/He ionization region Cepheid envelopes, con- 
vection reduces its source Y by an order of magnitude, but it is not capable of 
reducing it at the boundaries of the convective region to less than a value of 
about 2-3. 

3.4. Sequence of Cepheid models 

In Fig. |6] we display the behavior of the relative growth-rates r/ of the funda- 
mental and first overtone modes for a sequence of Cepheid models (M=5M Q , 
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Figure 5. Work integrands < pdv > for a typical Cepheid model in the 
middle of the instability strip; top: linear; and middle: nonlinear; 
bottom: dimcnsionlcss entropy gradient Y , the turbulent energy, normalized 
to unity, and the fraction of energy carried by the convective flux F c . 
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Figure 6. Behavior of the growth-rates (rj = 2k x Period) of the funda- 
mental and of the first overtone modes as a function of T e jj for a sequence of 
Cepheid models (M=5M©, L=1664L , X=0.70, Z=0.02) 

L=1664L , X=0.70, Z=0M). The r/'s are denned as rj = Attk/uj (where the 
linear eigenvalues are a = iu + n) and represent the pulsational energy gain per 
cycle (the inverse of the quality factor Q in electronic circuits). (The models 
for the four sequences differ slightly because the small, but nonzero turbulent 
flux affects the equilibrium.) The four recipes are seen to give almost the same 
results, again showing that the time-dependence of the recipes is not critical. 
This is also in agreement with what YKB found, namely that the sudden ap- 
proximation (instant adjustment of et to the instantaneous static solution) is an 
excellent approximation. This is in contrast to the frequently made, but bad 
approximation of 'freezing' the convective flux in the computation of the linear 
eigenvalues. 

We conclude this section by noting that the three recipes give very similar 
behavior, and that the apparent differences in the published results are most 
likely due to different choices of the a parameters than to the recipes themselves. 

4. Results 

4.1. Light-Curve Fourier Decomposition Coefficients 

In the Introduction we mentioned the presence of resonances in Cepheids be- 
tween the self-excited pulsational mode and an overtone, and that this is the 
reason why the Fourier decomposition coefficients of the light-curves and of the 
radial velocity curves show so much structure in the Cepheids. RR Lyrae are 
devoid of such resonances and show a relatively dull behavior of the Fourier de- 
composition coefficients as a function of period. The presence of these resonances 
puts very severe constraints on Cepheid models that any turbulent convective 
description must satisfy (see below). 




11 



Purely radiative models were unable to reproduce the large excursion of the 
"Z" shape of the fai Fourier coefficient for the overtone Cepheids (e.g. Buchler 
1998, Antonello & Aikawa 1993, 1995; Schaller & Buchler, 1994, unpublished 
preprint). It came therefore as a pleasant surprise when the nonlinear turbulent 
convective overtone Cepheid models showed the ability to reproduce wide excur- 
sions in the 021 such as indicated by the observations. However, since extensive 
nonlinear calculations are needed, a systematic study of the Fourier coefficients 
remains to be made. 

4.2. Double-Mode Pulsations in Cepheids and RR Lyrae 

The numerical modelling of steady, nonlinear double mode (DM) pulsations has 
been a long-standing quest in which purely radiative models have failed. In a 
recent paper (Kollath et al. 1998) it was found that with the inclusion of a ID 
turbulent convection recipe in the code, DM pulsations appear quite naturally in 
Cepheid models. Almost concomitantly, but independently, Feuchtinger (1998b) 
encountered DM pulsations in RR Lyrae, which we have since also confirmed. 

Buchler et al. (1999) show that the behavior of the double mode phe- 
nomenon both in Cepheids and in RR Lyrae can be captured by rather simple 
amplitude equations (also called normal forms). Here Aq and A± denote the 
modal amplitudes of the two excited vibrational modes. 

" ' l " (k - q Q0 A 2 - q 01 Al - TqAq) A (17) 



dt 

fMj («i-gio^-giiA? -r 1 Aj)A 1 (18) 



dt 

The transient evolution of a hydrodynamics model is governed by these equations 
(Kollath et al. 1999). Furthermore, and more importantly, the fixed points of 
these amplitude equations (dA^/dt = dA\/dt = 0) allow an overview of the 
'modal selection' (bifurcation diagram) in the physical space of L and T e ft, for 
example. 

Kollath et al. 1999 have computed double-mode pulsations with both the 
YKB and the GW formulations (Eqs. |l^). Both recipes give rise to DM pulsa- 
tions in broad domains of the physical (L, M,T e ff) and (a) parameter space. 
We note that the DM behavior occurs in Galactic as well as in MC Cepheid 
models (Kollath et al. 1999). 

We have recently found that correcting the recipe (Eq. |3| for small Peclet 
number (see below) allows one to find also overtone double-mode pulsations (as 
are observed), in which the first and second overtones are excited. 

We finish this section by stressing that the range of observed double-mode 
periods, amplitude ratios and Fourier phases provides a powerful set of con- 
straints on the model parameters. This is true not only for the DM fundamental- 
ly overtone and \ st -2 nd pulsators, but also for the single mode Cepheids and 
RR Lyrae. 

4.3. RR Lyrae pulsations 

Extensive computations of convective RR Lyrae model pulsations have been 
made by Feuchtinger (1998a, 1999) and Feuchtinger &: Dorfi (1997) using the 
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Gehmeyr- Winkler recipe (Eqs. [XOj ) , but omitting the turbulent flux and the 
turbulent pressure, and by Bono et al. (1997), using the Stellingwerf prescription. 
These results are reviewed by Feuchtinger in this Volume. 

5. Observational Constraints 

In the Introduction we have pointed out that purely radiative models have failed 
in many respects. In the previous section we have discussed the recent successes 
of the turbulent convective computations, and we have compared them to the 
results of the prior radiative modelling. Despite the early optimism that these 
successes have generated, there remain some very serious discrepancies. That 
they are more apparent in the Cepheid models than in RR Lyrae is not aston- 
ishing. RR Lyrae stars, despite their observed differences, are really very similar 
to each other in mass, luminosity and composition. The observational data of 
the Cepheids, on the other hand, span more than an order of magnitude in lu- 
minosity and in mass, and the metallicity of the Galactic Cepheids is four to 
five times that of the SMC Cepheids. These observations impose a number of 
important constraints that we now address. 

5.1. Resonances and Mass Luminosity relation 

In 1928 Hertzsprung noticed that a bump on the light curves of the funda- 
mental Cepheids that moved from the descending branch of the light curve to 
the ascending one at a pulsational period Pq ~ lOd. Later it was found that 
this Hertzsprung bump progression and its corresponding manifestation in the 
Fourier decomposition coefficients is related to a resonance of the fundamental 
mode of pulsation with the second overtone (Simon & Schmidt 1976, Buchler & 
Goupil 1984). 

It is well established by now that structure in the behavior of the Fourier 
decomposition coefficients of light curves and radial velocity curves as a function 
of, usually, the period or effective temperature, is related to resonances. Because 
the structure of the stellar envelope changes with this 'control' parameter, the 
periods and period ratios also change and the excited pulsation mode can run 
into a resonance condition with an overtone which then gets entrained by this 
resonance through nonlinear coupling (e.g. Buchler 1993). There are two reasons 
for the prominence of the Hertzsprung resonance. First, the period ratio P0/P2 = 
2 is small, making the modal coupling very low order and thus strong. Second, 
the second overtone, while linearly stable, is only slightly stable, which allows it 
to be entrained to an appreciable amplitude through the resonance. 

For the first overtone Cepheids, there is also a 2:1 resonance which occurs 
in the lightcurve data at a period around 3-4 d (Antonello et al. 1990). More 
recently this resonance has also been found to be quite prominent in the radial 
velocity data (Kienzle et al. 1999). 

The data from the microlensing surveys indicate that both these resonances 
occur at the same pulsation periods, to within a day. These constraints allow 
us to compute 'resonance' masses and luminosities as a function of Z for stars 
in the instability strip. In other words it allows us to determine points on the 
mass-luminosity curve that these stars must obey. In a recent Letter Buchler, 
Kollath, Goupil & Beaulieu (1996) discussed the implications of these resonances 
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for Cepheid models, and they showed that there is a serious discrepancy between 
the pulsational models and stellar evolution calculations, with in particular the 
pulsation masses being much too small for the SMC models. Those results were 
obtained with purely radiative models. 

We have since reexamined the problem with the turbulent convective mod- 
els. The computed resonance luminosities still turn out largely independent 
of Z as suggested by the observations (Beaulieu, private communication), and 
the resonance masses are still much smaller for the lower Z values. One notes 
though that lower masses are in agreement with the evolutionary calculations 
(e.g. Chiosi et al 1993, Schaller et al. 1992, Baraffe et al. 1998). However, rec- 
onciling these results with reasonable widths for the instability strips remains a 
problem. 

5.2. Largest First Overtone Period 

The largest observed periods for the Galactic overtone Cepheids are Pi max = 
7?57 for V440 Per, P lmax = 5?44 for X Lac (Antonello et al. 1990). Beaulieu 
quotes Pi max = 5<?84 for star No. 114 in the LMC, Pi max = 5^98 for star No. 255 
in the SMC. The star V440 may be an oddball and if we can ignore it, then there 
appears to be a metallicity independent upper limit Pi max ^ 6d. The overtone 
periods in the beat Cepheids do not provide any upper limit here, because all 
of them are considerably smaller than the periods of the single-mode overtone 
Cepheids. 

For the following discussion we refer to Fig. 12 of YKB that shows Hertz- 
sprung -Russell diagrams with the shapes of the fundamental and first overtone 
instability strips (IS) for Cepheids. The overtone IS is pinched off above a cer- 
tain luminosity (or period). We refer to the overtone period at the tip of the 
overtone IS as Pi max . 

We can adjust the a parameters so that for Galactic Cepheid models (typi- 
fied by X=0.70, Z=0.02) we get an upper limit of Pi max = 6d. If, with the same 
a's we now compute the P\ m ax f° r SMC composition (typified by X=0.726, 
Z=0.004) we obtain values exceeding 12 d for a broad range of a parameters, 
which is clearly not in agreement with the observations. Furthermore, if we 
take a larger Z value for the Galaxy, as is sometimes suggested, the discrepancy 
worsens. Conversely, adjusting the a's to get Pi max in agreement with the SMC 
value gives maximum overtone periods that are far too small for the Galaxy. 

The problem arises because the lower masses and the low metallicity of the 
MC Cepheids. For the same luminosity and a's these stellar models are more 
unstable (larg er t/'s) than their Galactic siblings, which pushes the Pi max to 
much higher values. At this time a resolution of this difficulty appears to pose 
a serious challenge to convective Cepheid modelling. 

5.3. Width of Instability Strip 

Adjusting our a's so that the Galactic Pi max is ~ G^O we invariably obtain (at 
the luminosity for that overtone period) a width for the fundamental IS that 
exceeds 1000K. We consider such a width quite excessive because in addition to 
the intrinsic width (for a given composition) there must be substantial widening 
due to metallicity dispersion. 
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We have found that this problem may have its origin in the formulation 
of the convective recipes Eqs. [l(], that are all derived under the assumption of 
efficient convection, i.e. of large Peclet number (Pe is the ratio of radiative to 
convective diffusion times, thus measures the radiative damping of the convective 
elements). But this assumption of large Pe is not satisfied everywhere in the 
Cepheid envelope. In Fig. we display the Peclet number throughout typical 
Cepheid envelopes, namely M=5M Q , L=5OOOL , with T e yj=4600K (solid line) 
and T e j-y=5300K (dotted line). As a guide we also show the behavior of the 
dimensionless entropy gradient Y in the upper figure. The Peclet number is 
very large throughout most of the combined H/He 1 partial ionization region, 
but as the blown up bottom picture shows, the second helium ionization region 
lies entirely in the low Pe regime. There is also a small surface region, below 
T=8000 K, where Pe < 1. Since the second ionization stage of helium contributes 
substantially to the work-integrand we expect that an improved treatment of 
ineffective convection has a major effect on the growth-rates. 

In order to accommodate the small Pe limit CD suggest an interpolation (see 
also KuhfuB 1986), which after a small manipulation can be seen to be equivalent 
in the GW scheme to limiting J, i.e. both F c and S in Eqs. ( JlOj ) by a factor of 
the form 

f - I 

J** ~ i + ar p e -i ' 

where the Peclet number is the ratio of convective and radiative the diffusion 
coefficients, Pe= D c /D r , with 

_ 4 acT 3 

and 

D c = AeJ /2 

The definition of Pe is somewhat arbitrary and we have introduced an addi- 
tional free, order unity, parameter a r . Both the convective flux and the source 
term thus become proportional to Pe when Pe is small, and scale with an addi- 
tional factor of e\ /2 . As a result, in the static equilibrium model (if we disregard 
diffusion), one obtains et ~ Y 2 , instead of ~ Y, and both F c and S scale with 
Y 3 , instead of Y 3/2 . 

Referring back to Fig. [5], we note that in the higher mass Cepheid models, 
such as this one, a strong (Y very small) convective zone can appear in the Fe 
ionization region ~200,000 K. 

We have performed some preliminary calculations with the Peclet factor 
fpec- Figure || shows the change of some of the model properties for M=5M Q , 
L=3184, A=0.70, Z=0.02, with a r =0, 0.1 and 0.25. (a c = 4.0, a v = 1.0, a s 
= 0.75, a p = 0.667, at = 1.0, ad = 4.0, «a = 0.58). One sees the reduced effi- 
ciency of turbulent energy source particularly in the H/He region. This figure is 
somewhat misleading though. The Peclet factor decreases the turbulent energy 
which means that in order to satisfy the observational constraints, we have to 
compensate by increasing some of the a's, which is not reflected in the figure. 
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5.5 5 4.5 4 

Log T 

Figure 7. Top: Entropy gradient (Y); Middle: Peclet number, as a func- 
tion of Log T for Cepheid models: M=5M , L=2OOOL , X=0.70, Z=0.02; 
T e// =4900K (solid line), T e// =5300K (dotted line); Bottom: Blow-up of 
small Fte number scale. 

However, the important result is that the introduction of f pec has a differ- 
ential effect on the stability of the fundamental and first overtone modes that 
goes in the direction of solving the problem of the excessive width of the fun- 
damental instability strip. That width is most sensitive to the parameter a r 
which can thus quite satisfactorily be used for adjusting the width to satisfy 
this observational constraint. 

5.4. RR Lyrae 

For RR Lyrae the 'linear' constraints are somewhat different. Here, because glob- 
ular clusters have approximately the same luminosity, we know the edges of the 
corresponding instability strips directly from the observed fundamental (RRab) 
and first overtone (RRc) periods. The observed T e ff are probably too uncertain 
to be directly very useful, although for specific globular clusters stars they may 
provide useful constraints on the temperature widths of the ISs. However, it 
is important to keep in mind that using the linear periods for the constraints 
can be misleading because a model may be inside the linear instability strip, yet 
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Figure 8. Effect of Peclet correction factor on, from left to right, top 
down, as a function of zone index, the entropy gradient Y, the turbulent 
energy e*, the convective flux (F c /Ftot), the temperature versus zone index, 
the turbulent flux F t /F to t, the work- integrands of the fundamental and first 
overtone modes. 

Solid-dotted line: a r =0, dotted lines: a r —0.1 and thin line: a r ~0.25. 



the corresponding limit cycle may be unstable (e.g. Buchler & Kovacs 1986). 
In other words, the observational instability strip is narrower than the linear 
instability strip, and it takes nonlinear calculations to determine the former. 

A great deal of observational lightcurve data are available, and their Fourier 
decomposition coefficients provide strong nonlinear constraints that need to be 
satisfied by the models. In particular there is a nice correlation between the 
coefficients that has been discussed by Kovacs and Kanbur (1997). 

The double-mode, RRd stars also provide tight constraints on the alpha 
parameters through the observed periods and period ratios, and through the 
amplitude ratios. 

5.5. Additional Constraints 

From the observational data one can infer period - radius relations for Cepheids 
(e.g. Laney & Stobie 1994). These constraints are very basic since they con- 
cern the equilibrium models. Bono et al. (1998, 1999) have had mixed success 
reproducing this observational constraint. 

Other nonlinear interesting constraints arise from the temperature fluctu- 
ations of the Cepheids and RR Lyrae. Simon, Kanbur & Mihalas (1993) have 
shown that essentially all Cepheids have the same T max (same spectral type), 
but that T m i n depends on the amplitude of pulsation. Similarly, but oppositely, 
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RR Lyrae have essentially the same T m i n (e.g. Kanbur & Phillips 1996). These 
constraints have not yet been taken into account in the calibration of the a 
parameters. 

6. Theoretical Difficulties 

For a given set of as the ID recipes all give the same turbulent energy et in the 
equilibrium models when the turbulent flux is neglected, which is 

3T 

e t = a 2 s a\ — c 2 s Y or u t ~ c s W (19) 
J- 1 

When this et is inserted in the expression for convective flux one obtains 

( pCp T)c s Y 3 / 2 (20) 

However, F c is physically limited to transporting the available energy fluctua- 
tions whose very generous upper limit is (pc p T). Furthermore, under the under- 
lying assumptions the convective velocity should be subsonic. Eq. (|20|) clearly 
exceeds this upper limit when Y is larger that unity (cf. Fig. ^). For this reason 
ad hoc flux limiters have been proposed (e.g. Feuchtinger 1999) 

The ultimate reason for the breakdown of the ID recipe is twofold. First, 
in the convective boundary regions where Y ^ 1, the turbulent Mach number 
Mt = (ut/c s ) invariably becomes important and can even exceed unity in some 
weakly convective models. However, it is well known that the approximations 
leading to the recipes Eqs. (g) break down, in particular because pressure fluc- 
tuations were neglected in their derivations and the dissipation term e does not 
take into account larger Mach numbers. Second, as stressed by Canuto (1998) 
the downgradient approximations for the fluxes (here F c ~ ds/dr) are not very 
satisfactory. Neither of these problems can easily be corrected within a practi- 
cally useful ID recipe, and a flux limiter may be the most expedient available 
patch. 

7. Conclusions 

The introduction of a simple recipe for turbulent convection, namely a ID time- 
dependent diffusion equation for the turbulent energy and the concomitant ex- 
pressions for convective and turbulent fluxes, eddy viscous pressure and turbu- 
lent pressure, have provided a substantial improvement over purely radiative 
models. Perhaps the most remarkable achievement is the successful modelling 
of double-mode pulsations both in Cepheids and in RR Lyrae. Correcting the 
recipes for low Peclet number provides a further definite improvement. 

However, a number of problems persist as we have seen. The number (8) 
of free, order unity parameters is large which makes a search for optimal values 
quite burdensome. We still hope that it will be possible to find a range of suitable 
parameters that satisfy all the observational constraints, and this independently 
of metallicity or of the stellar type. If not we will be forced to adopt more 
complicated formulations of convective transport. 
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